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Abstract 

Duality is considered for the perturbation theory by deriving, given a series 
solution in a small parameter, its dual series with the development parameter 
being the inverse of the other. A dual symmetry in perturbation theory is 
identified. It is then shown that the dual to the Dyson series in quantum 
mechanics is given by a recently devised series having the adiabatic approx- 
imation as leading order. A simple application of this result is given by 
rederiving a theorem for strongly perturbed quantum systems. 
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A known result of fluid mechanics jl], given the Navier-Stokes equation 
in absence of any forcing 

— + (u • V)u - i/V 2 u + Vp = 
at 

where v is the viscosity, p the pressure and u the velocity field, is that, by 
fixing the 0-th order solution through the Eulerian part as [|| 

du 



dt 



+ (u • V)u = 0, 



one obtains a perturbation series for a large Reynolds number Re while, taking 
at the leading order the equation 

— — - z^V z w = 
at 

one obtains a perturbation series for small Re. So far, this was the only case 
in perturbation theory, as applied to physics, with an equation generating 
both a small perturbation series and its strong perturbation counterpart. 

Duality in perturbation theory should then be understood in the sense 
that, for a given differential equation, one has the possibility to derive pertur- 
bation series both in a small development parameter and its inverse, giving 
in this way the possibility to study a solution in different regions of the pa- 
rameter space. However, it should be said that the method could not be of 
absolutely general usefulness as some limitations can appear for the compu- 
tation both of the leading and higher orders. Beside, there exists situations 
where better approximations are known, as I will show. However, it is easy 
to realize that a lot of problems in physics can get new insights from this 
approach, so that, it is worthwhile to be exploited. 

A natural question, in the light of the above defined duality in perturba- 
tion theory, is what should be the dual to the well-known Dyson series for the 
Schrodinger equation. The answer to this question is the main aim of this 
communication. In fact, the existence of this possibility gives a new technique 
to analyse quantum systems in different regions of their parameter space. 

Quite recently, I showed that new solutions for the Schrodinger equation, 
in time-dependent problems, can be obtained when a strong perturbation is 
applied to a quantum system Q. This approach seems to indicate that, at the 
leading order, an adiabatic approximation should be used J3j. On a different 
line of research, Mostafazadeh || was able to show that a series exists, for 
the Schrodinger equation, with a well defined development parameter, having 
the adiabatic approximation as leading order. Using duality in perturbation 
theory the above different research lines can be merged, as I am going to 
show, giving us the main result of this paper. In fact, the series derived 
by Mostafazadeh is dual to the Dyson series. Then, the theory of strong 
perturbations in quantum mechanics can be proved to be dual to the standard 
small perturbation theory. 

In order to show how a dual series can be obtained in a simple case, let 
us consider the model given by the following differential equation 



x = fo(x) + Xfi(x) 



where the dots mean derivation with respect to the time and A is an ordering 
parameter. It is a well-known matter that, when A — > 0, a solution series 
of the form x ~ xo + Axi + A 2 X2 + 0(A 3 ) can be obtained. However, as for 
the Navier-Stokes equations, we are free to choose at the leading order, as an 
unperturbed equation, xq = A/i(xo). To show that this choice gives a dual 
perturbation series, I rescale the time variable in eq. (|) as t -> V\t = t. One 
gets 

Ax = f (x) + A/i(x) 

where now the dots mean derivation with respect to r. It is quite easy to 
verify that the series 

x = x + t x i + To x 2 + O 
A X z 

is a solution of eq. (0) when 

x'o = fi(x ) 

x'i = f[(x )x 1 + /o(x ) 
x'2 = fo{x )xi + f[{x )x 2 + ^/a (x )xj 




By analogy with the results in quantum mechanics |J, I take the above as 
the dual method to small perturbation theory to obtain a dual perturbation 
series to a given one. It is important to note that the above result is true 
independently of one's ability to solve the leading order equations. 

We see that the arbitrariness in the choice of the leading order equation 
gives rise to a symmetry. In fact, putting A = 1 into eq.(|j), there is no more 
reason to see any difference between the perturbation and the unperturbed 
system in the same way as happens in fluid mechanics. This means that the 
series given by the small perturbation theory can be derived from the one given 
by the dual method and vice versa by simply interchanging fo(x) <-> /i(x). 
That is a symmetry of the perturbation theory whose meaning can be really 
understood only after the introduction of the dual series. Actually, the general 
solution of eq.(||) can be written, for the one-dimensional case, as 

f x i 1 
t — to = / dx' = 

J *o V2^E + f* Hx")dx" + h{x")dx" 

with E a motion constant. It is easily seen that both the series expansions, 
for small /o or f\, can be straightforwardly obtained. What it is interesting is 
that, the small parameter in a case is the inverse of the development parameter 
in the other. From the discussion above it should be clear that both series 
can have the same problems as secularities or divergent terms. 



A more interesting example is given by the following Duffing equation 



X + UJqX + f3x i = fo cos(idt). 



■ 2 0fi 



By setting r = uJot, v = £, = -j^x and A = -^§-, one gets the rescaled 
equation 

| + £ + Xf = cos(zat). 

where the dots mean derivation with respect to r and A is just a parameter 
measuring the strength of the nonlinearity. Eq. (|l0|) is generally considered, 
analytically, only for small A, but what happens in the limit of a very strong 
nonlinearity? Duality can be applied and one easily realizes that, for large 
values of the parameter A, the quantity e = + A^£ 4 tends to be a constant 
of motion. This is due to the result that the perturbation completely drives 
the system. So, we have regular periodic motion in the considered limit. This 
example shows that, although the leading order equation can be easy to solve, 
going to higher orders could be very involved. 

A class of important problems arises from the Schrodinger equation that 
I consider in the one-dimensional form 

where A —* oo. One could apply immediately the symmetry between the dual 
and small perturbation theories discussed so far and use without difficulty 
the Rayleigh-Schrodinger approximation scheme. While that is a correct ap- 
proach, I will show how the dual method works in this case. So, let us put 
£ = y/Xx, ip = + + O (^) and E = \E + E x + O U). This yields 



the following equations 

h 2 c^Vo 



+ Fi(e^o = Wo 

- + ( f p 

2m d£ 



h 2 d 2 ^ 



2m d£ 2 

- + Vi(eOV>i + Vo*Po = ExipQ + E i>i 



where e = At the leading order we get a well-known equation, that is, a 
second order differential equation with a slowly varying coefficient due to the 
perturbation. In this case we can apply the WKB approximation J?]]. Thus, 
the dual method yields in this solution that is a combination of both 

Rayleigh-Schrodinger and semiclassical methods. 

There are several problems where the above approximation can be applied. 
A well-known example is given by the anharmonic oscillator that has a large 
body of literature || and is a model which any approximation scheme should 
address. The Hamiltonian can be cast in the form 



The method I discussed so far gives an unambigous answer to this problem, 
i.e. the leading order approximation, when the anharmonicity is very strong, 
can be obtained by solving the equation 



\^ + \ q ) Mq) = eMq)- 

The quartic oscillator is well-known in the literature ||. So, we can compare 
our method with numerical results. To leading order of the WKB approxima- 
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tion of the energy levels, normalized in unit of \ \-\J\ ^ j the agreement is 

within 18% with the true value of the ground state energy for the anharmonic 
oscillator. That agreement improves for higher excited states. However, we 
know from Symanzick scaling that the quartic oscillator is the right approxi- 
mation for energy levels of the anharmonic oscillator when A — > oo [jlCfl . Often, 
the use of semiclassical eigenfunctions can be too much involved and better 
approximation schemes, as those given in ref. ||, can improve the situation. 
Duality, as applied in perturbation theory, yields anyway a definite answer. 

The situation is surely more interesting in time-dependent problems. By 
noting that the only meaningful quantities are transition probabilities between 
states of the unperturbed system, we have the initial conditions definitely fixed 
breaking in this way the dual symmetry of the perturbation theory. In fact, 
in ref. |J I showed that the problem 

(H + \V(t))\iP)=ihd t \i;) 

d 

with dt = — and A — > oo, using the above dual method, has the leading order 
at 

solution 

U{t)\n) + 0' 1 



with 



XV(t)U(t) = ihd t U(t) 

and Hq\u) = E n \n). So, the unperturbed solution fixes the initial condition 
as also happens in the small perturbation theory. But, in order to leave the 
dual simmetry untouched, one should physically consider also systems ini- 
tially prepared with the eigenstates of the perturbation, but this is not the 
case for the computation of probability transitions. Then, we can conclude 
that, for the time-dependent perturbation theory as usually applied in quan- 
tum mechanics, the dual symmetry is broken due to the choice of the initial 
conditions. 

By the methods discussed above, we can obtain the main result of the 
paper. Our aim is to show that the dual to the Dyson perturbation series, for 
the time-dependent Schrodinger equation, is given by the series obtained by 
Mostafazadeh S having the adiabatic approximation as leading order. So, 



let us consider the Schrodinger equation H(t)U = ihdfU being U the time 
evolution operator. The Dyson series is the solution of that equation and can 
be written formally as 

U(t) = Texp (-jrf* dt'H(t'^j 

with T the time-ordering operator. This is a compact writing for the series 
development 

17(f) =1-1 f dt'Hit') + (- 1 -) 2 C dt' f dt"H(t')H(t") +••• 
n Jo \ h) Jo Jo 

A dual series to the one above is meant as a series having for develop- 
ment parameter its inverse, as discussed above. To obtain it, let us con- 
sider the case with the hamiltonian H being constant in time. Assuming, 
for the sake of simplicity, here and in the following that the hamiltonian 
has a discrete spectrum, the solution to the time-dependent Schrodinger 
equation is easily obtained through the time evolution operator U(t) = 
^2 n e~K Ent \n)(n\ with H\n) = E n \n), where the effect of the time deriva- 
tive is simply dfU(t) = J2 n ( — e~~a Ent \n) (n|. Instead, for a time- 
dependent hamiltonian H(t), in the case of the adiabatic approximation we 
have U A (t) = J2 n e ian(t) \n,t){n,0\ being a n (t) = 7n (f) - ^^dt'E n {t') with 
H(t)\n,t) = E n (t)\n,t), so that E n (t) gives the dynamical part of the phase 
a n (t) and j n (t) = i J dt' (n,t'\d t '\n,t') the geometrical part. It is natural to 
ask how one can define a derivative Dt to obtain the same result as for the 
time-independent case, that is, D t UA(t) = J2 n e ian ^\n,t)(n,0\. 
It is quite easy to verify that the following definition of Dt has the requested 
property 

Dt = dt + i ^ (m, t\idt\n, t)\m, t)(n, t\ 

so that, the adiabatic approximation is exact for the equation 

H(t)U A (t) = ihD t U A (t) 

and it is easily verified that iDt\k,t) = 7fc(f)|A;,f}. In this way, we are a step 
away from the sought result. In fact, let us introduce a generic perturbation 
V(t) into eq.©, so that 

(H(t) + V(t))U(t)=ihD t U(t). 

In this form the duality principle of perturbation theory can be ap- 
plied. Having V(t) as the unperturbed hamiltonian and H(t) as the per- 
turbation, we get no physical meaningful results from the leading order 
equation V(t)U^(t) = ihD t U^\t) unless we choose V(t) = H(t) or 
V(t) = — (m,t\ihdt\n,t)\m,t)(n,t\. The latter is the interesting case giv- 

n, m 

ing trivially the standard Dyson series. So, when we take H{t) as the un- 
perturbed hamiltonian and V(t) as the perturbation, we have at the leading 
order 
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H(t)U^(t) = ihD t U^(t), 



but this is nothing else than eq,(|2~ll) and then U^°'(t) = U^it) (i-e. at the 
leading order we have the adiabatic approximation) . To complete the identifi- 
cation with the Mostafazadeh result, we have that the higher order corrections 
are computed solving the equation 

H'(t)U'{t) = ihd t U'(t), 

where 

H'{t) = U ] A {t)V{t)U A {t) = - e- i{amit) ~ an{t)) {m,t\ihd t \n,t)\m,0){n,0\ 

n,m,n^m 

as given in ref. ||. Then, one obtains ||] 

U(t) = U A (t)Texp (-^^dt'H'it'^j 

that completes the proof. No adiabatic hypothesis entered into this argument 
as it should be. 

As an application of that result, a theorem recently derived by me and, 
in a rigourous way but in a different context, by Joye Q can be obtained for 
the theory of the strong perturbations in quantum mechanics ||. In fact, it 
was proved that, for a quantum system described by the Schrodinger equation 
(Hq + XV(t))\ip) = ihdt\ip), in the limit A — > oo, the adiabatic approximation, 
using the eigenstates of the perturbation V(t), is a good approximation for 
\ip). So, by considering the perturbed system in the interaction picture, gives 
the hamiltonian Hi(t) = U (0 ^ (t)XV(t)U^ (t) where H U<®(t) = ifid t U^(t). 
Then, the result obtained for the Dyson series by duality can be directly 
applied to Hj{t). We obtain for a small A (otherwise we miss convergence), 
the Dyson series 

U I (t) = I- % -J^dt'H I (t') + ( K -^j 2 J^dt'J^ dt"H I (t')H I (t") + --- 

and, for a large A, the Mostafazadeh result applied to Hj(t) having at the 
leading order the adiabatic approximation as it should be. It must be noticed 
that, in the latter case, the eigenstates to be considered are those of the 
perturbation. In fact, we have C/(°) t (t)y(t)LT(°)(t)|A;,t)/ = vl(t)\k,t)i but 
this is equivalent to V(*)(l7W(t)|M)j) = ^(*)(^ (0) (*)|M)/)- Then, v[ is 
an eigenvalue of V{t) and U^°\t)\k, t)j differs just by a time-dependent phase 
factor from the corresponding eigenstate of the perturbation. As a by product, 
we get the confirmation that higher order corrections are those computed by 
the method of strong perturbation theory as given in Q. 

In summary I introduced the duality principle in perturbation theory for 
differential equations. A dual method with respect to the theory of small 
perturbation is yielded and a dual symmetry between the two methods arises 
from the freedom in the choice of what the perturbation is. The use of dual- 
ity shows that, for the time-dependent Schrodinger equation, the dual series 



to the Dyson one is given by a perturbation series computed recently by 
Mostafazadeh, with a well-defined development parameter, having the adia- 
batic approximation as leading order. This enriches the possibility to analyze 
quantum systems in completely different regions of their parameter space. 
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